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Abstract. We discuss the role of noise and dissipation in the explosive spin- 
odal decomposition scenario of hadron production during the chiral transition 
after a high-energy heavy ion collision. We use a Langevin description inspired 
by noncquilibrium field theory to perform real-time lattice simulations of the 
behavior of the chiral fields. Preliminary results for the interplay between ad- 
ditive and multiplicative noise terms, as well as for non-Markovian corrections, 
are also presented. 
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1. Introduction 

For high enough values of temperature, strongly interacting matter should be in a 
quark phase due to asymptotic freedom. In fact, finite-temperature lattice QCD 
simulations provide strong evidence of a deconfined quark-gluon plasma (QGP) 
phase at sufficiently high temperature Moreover, relativistic high-energy heavy 
ion collisions can probe strongly interacting matter under extreme conditions. They 
provide valuable information on the new state of matter that seems to have been 
created according to recent data from experiments at BNL-RHIC [EJ, even if its 
true nature is still uncertain. 

As the QGP presumably created in a heavy-ion collision expands, it cools down 
and, at T c ~ 150 MeV according to the Lattice [P, undergoes a phase transition 
(or a crossover) back to hadronic matter. Results from CERN-SPS and BNL-RHIC 
feature what has been called sudden hadronization [ 3 or explosive behavior [0] in 
the hadronization process of the expanding QGP and seem to favor a fast (explosive) 
spinodal decomposition scenario as the mechanism of phase conversion. This should 
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be confronted to the process of bubble nucleation, whose time scales are larger due 
to the existence of a barrier to overcome [ 5 . Possible signatures of the explosive 
behavior in high-energy nuclear collisions were proposed in Ref. 

Effective field theory models for the chiral and the deconfincment transitions 
in QCD also indicate that most of the plasma is quenched into the spinodal region, 
due to the fast expansion of the system as compared to the low nucleation rate, 
then undergoes an explosive phase conversion [0| El EH El This leads to 
what we will refer to as the explosive spinodal decomposition scenario. Although 
an expansion term in the evolution equation for the order parameter plays a role 
somewhat similar to dissipation, most studies do not take into account genuine 
dissipative effects due to the interactions with the medium. However, dissipation 
effects have proved to be important, for instance, in the context of disoriented chiral 
condensate (DCC) formation in heavy ion collisions [El El El El El Ej and 
could be an important ingredient in the explosive spinodal decomposition scenario. 

Recently, we have considered the effects of dissipation on hadron production 
during the QCD transition after a high-energy heavy ion collision in the simplest 
fashion [EJ- Using a phenomenological Langevin description for the time evolution 
of the order parameter, inspired by microscopic noncquilibrium field theory results 
[ El H41 1191 |2*U| , we performed real-time lattice simulations for the behavior of 
the inhomogeneous chiral field. We concluded that the effects of dissipation could 
be dramatic even for very conservative assumptions. It was shown that even if the 
system quickly reaches this unstable region there is still no guarantee that it will 
explode. 

Nevertheless, the description implemented in Ref. [El captures only a small 
part of the much richer spectrum of possibilities for dissipation and noise effects 
in the dynamics of phase transitions. Part of this richness is captured by a lin- 
ear response nonequilibrium quantum field theory analysis, which usually provides 
not only different contributions to noise and dissipation terms, but also compli- 
cated memory kernels instead of simple Markovian terms proportional to the time 
derivative of the field [El El El HOI 

In what follows, we present our findings from a traditional Langevin analysis of 
the evolution of the order parameter for the chiral transition, and discuss prelimi- 
nary results including multiplicative noise [ I21| and memory effects [E2- A robust 
feature that comes about from this study is that all the contributions tend to add 
up to delay the process of thermalization. 

2. Dissipation versus explosive decomposition 

We consider the real-time dynamics of chiral symmetry breakdown of a QGP created 
in a high-energy heavy ion collision [El- We assume the system to be characterized 
by a coarse-grained free energy 
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where V e f / {4>, T) is an effective potential of the Landau-Ginzburg form whose co- 
efficients depend on the temperature, and <fr(x, t) is a real scalar field which plays 
the role of an order parameter that is not conserved, such as the chiral condensate. 
To model the mechanism of chiral symmetry breaking found in QCD, we adopt 
the linear a-model coupled to quarks, whose standard Lagrangian can be found, 
for instance, in Ref. [ 110] . This approach is widely used in the literature and its 
specificities imply no major limitations to our main results. Quarks are integrated 
out yielding the effective potential V e ff(<fi,T). The pion directions play no major 
role in the process of phase conversion we are considering, so we focus on the sigma 
direction represented by the field cj) [10.. For simplicity, we ignore effects due to 
the finite size of the plasma [I23|. 

The framework for the dynamics is assumed to be given by the following 
Langevin equation: 

__ V ^ + r? _ + __ = ^, t ), (2 ) 

where ^ is a real scalar field and 77, which can be seen as a response coefficient that 
defines time scales for the system and encodes the intensity of dissipation, is usually 
taken to be a function of temperature only, rj = i](T). The function £(a?, t) represents 
a stochastic (noise) force, assumed Gaussian and white, so that (£(x, t)) = and 
(£(x,t)£(x',t')) = 2r)T5(x — x')S(t — t'), according to the fluctuation-dissipation 
theorem. 

In our numerical simulations we solve Eq. JSJ on a cubic space-like lattice with 
64 3 sites under periodic boundary conditions, with a lattice spacing of a = 0.91 fm. 
We use a semi-implicit finite-difference scheme for the time evolution and a Fast 
Fourier Transform for the spatial dependence [ |^] . Temperature is fixed to the 
spinodal value T sp as 108 MeV [E3]- We perform several runs starting from different 
random initial configurations around the inflexion point of V e f / which happens at 
4>o 0.162 T and then average the results from the different initial configurations. 
For time steps of At = 0.001/T the results become independent of the lattice 
spacing once it is smaller than a ~ 1 fm. 

We show results of simulations for three different values of the dissipation coef- 
ficient, namely r]/T = 0, 2 and 4. It can be argued that the response coefficient has 
the form rj(T) ~ 2T/b, where b is a number of order one to first approximation [ 
I25| . The cases considered provide a conservative band around the value r/(T) w 2 T 
to illustrate the effect of dissipation. 

In Figure^we show the average value of <j> in units of its vacuum value, <p vac , as 
a function of time for the three different values of r\ mentioned above. The results 
clearly show that even for a very conservative value of dissipation, 77 = 2 T, the 
effect can be dramatic. For this value of 77, dissipation retards the time evolution 
of 4> towards its vacuum value in ~ 100% compared to the case with r\ — 0. The 
important point to be noted here is that for expansion times of the order of 5 fm/c, 
which is of the order of the time scales for RHIC collisions, there might be not 
enough time for the onset of the spinodal explosion. 
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Fig. 1. Average value of the chiral field tfi in units of its vacuum value cf> V ac as a 
function of time for rj/T = 0, 2, 4. 

3. Improved Langevin approach 

Although phenomenological equations of the form of Eq. (|2J) are widely used in the 
literature, formal derivations of effective equations of motion based on a more com- 
plete field-theoretic description of nonequilibrium dissipative dynamics EUl 120] 
show that much more complicated equations emerge. In general they are non-local 
equations of motion, with colored, non-Markovian dissipative and noise kernels, ex- 
hibiting both noise and dissipation terms that depend on the field amplitude. It 
is then important to analyze the possible outcome of dynamics emerging from the 
more realistic equations as compared to the case of the simpler phenomenological 
equation @. In fact, on very general physical grounds, one expects that dissipation 
effects should depend on the local density ~ (f) 2 <j> and, accordingly, the noise term 
should contain a multiplicative piece ~ <f>. These results emerge as approximations 
to the effective equations of motion derived e.g. in Refs. [O ED EH- 

Motivated by the results referred above, we consider an improved Langevin 
framework which includes the effects of multiplicative noise and density-dependent 
dissipation terms in the description of the time evolution of the order parameter. 
To assess the effect of the new contributions on the relevant time scales for phase 
ordering, one can adopt a scalar X(j> theory in the broken phase, whose standard 
Lagrangian density has the form: 



(3) 
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where </> is a scalar field playing the role of the order parameter and A > 0. We 
perform (3+l)-dimcnsional real-time lattice simulations to study the behavior of the 
inhomogeneous scalar field, taking into account lattice counterterms that guarantee 
lattice-size independence [El EH EH- 

To define an initial value problem for the evolution of the order parameter, one 
can assume that the system was quenched from a high-temperature environment 
in which the sign of the mass term was the opposite, corresponding to a symmet- 
ric phase, down to a temperature below some critical temperature of symmetry 
restoration (T < T c ), in which the effective potential has the form shown in the 
Lagrangian density [3] 

In our analysis, the time evolution of the field <j)(x, t) at each point in space and 
its approach to equilibrium will be dictated by an improved Langevin equation of 
the form 



^ - V 2 <j>(x, t) + [ m + m 4> 2 (x, t)] ' + V'{4>) = & {x, t) + cj>(x, t)& (x, t) , 



which models the local approximation for the effective equation of motion derived 
e.g. in Ref. [Ej. In Eq. r\\ and 772 will be taken to be functions of temperature 
only, rji = T)i(T). The functions £1(2;, t) and £2(2^ i) represent stochastic (noise) 
forces, assumed Gaussian and white. 

Eq. Q could, in principle, be obtained from a microscopic field-theoretic de- 
scription of the real-time nonequilibrium dynamics of the chiral field at finite tem- 
perature [ 1141 119j . The noise and dissipation terms, which originate from quantum 
fluctuations, are engendered by either self-interactions of the chiral field or coupling 
to one or more different fields that play the role of a heat bath, provided one incor- 
porates higher-order terms in the computation of the effective equation of motion 
for <p(x,t). In fact, it is well-known that one has to go up to two- loop corrections 
in order to pick up imaginary parts in the self-energy associated with viscosity and 
dissipation [ 1141 llfil 1191 I2()| . Self- interactions of the 4> field, as well as possible in- 
teractions with other fields in the medium, fully justify the inclusion of dissipation 
and noise terms such as done in the framework adopted here. 

In our numerical simulations we solve Eq. Q on a cubic space-like lattice. We 
use the leap frog method (see e.g. Refs. [|2U E3) f° r the time evolution and treat 
the Laplacian using a Fast Fourier Transform for the spatial dependence [E3|- To 
minimize lattice artifacts, we use counterterms calculated in lattice perturbation in 
the classical theory (tadpole plus setting-sun diagrams) [ I28j : 



where T is the temperature, a is the lattice spacing and /i is an arbitrary renormal- 
ization scale. It should be noted that the addition of this counterterm is supposed 
to guarantee independence of the results on a at equilibrium only. We checked 
the lattice independence of the results [ I21| , and adopt units such that m = 1 in 
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Fig. 2. Average value of the field cf> as a function of time under the effect of additive 
and multiplicative noises — 1 and 772 = 1, 2, 4) 

the Lagrangian of Eq. J2J), i.e. all dimensional quantities are expressed in units of 
m. In our simulations, the lattice extent is fixed to L = 16, and the time step is 
At = 1CP 3 ; T = 1, A = 1, and fj, = 1. The initial configuration is taken to be a ran- 
dom distribution of the form <j>(x, t) — (j> sp + 0.001 * rand(seed), where <f> sp = l/v3 
is the spinodal point of the bare double- well potential, and rand(seed) is a uniform 
random number distribution in the interval [—1,1]. At the end, results are averaged 
over several runs for different realizations of the noises. 

In Figure we present preliminary results for the volume average of as a 
function of time for different values of the multiplicative dissipation parameter 772 • 
For simplicity, 771 and 772 are taken to be constants. More realistic paramctrizations 
will appear soon [EI]- As 772 increases, the effect is of delaying the equilibration of 
the order parameter, as expected. 

4. Non-Markovian corrections 

As discussed before, the structure of memory integrals and colored noise that ap- 
pear in a realistic field-theoretic description of the dynamics of phase transitions is 
often rather complicated [EH1 00]' To develop and test new approaches, one can 
consider the much simpler case of dissipation in quantum mechanics, where all the 
approximations and important scales are under control, and where one also finds 
inumerous applications [I31|. 
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Starting from the noncquilibrium evolution of a particle coupled linearly to 
a set of harmonic oscillators in the Caldeira-Leggett fashion [ |2] , we study the 
effects of the non-local dissipation kernel, as well as of colored noise, that appear in 
the complete Langevin equation for the particle coordinate in space [ I22|. There, 
the memory kernel has its origin in the Feynman influence functional of the heat 
bath. To approach the kernel in a simpler, analytic way, we develop a systematic 
expansion in time derivatives whose convergence is regulated by increasing powers 
of the frequency cutoff in the distribution of oscillators, f2. Using this method, one 
can incorporate consistently corrections coming from the memory integral. The 
equation of motion thus obtained has the following form [I22|: 



which reduces to the traditional Langevin equation with white noise in the limit 
ilt — > oo, since I n (oo) — * 5 n oTr/2 and the noise correlator tends to 2rfT8(t — t'), 
consistently with the fluctuation-dissipation theorem. /„ (A) are integral coefficients 
which can be expressed in terms of incomplete gamma functions. Inspection of 
(|B"|) shows that terms containing higher-order time derivatives of Q are strongly 
suppressed by increasing powers of l/f2. 

5. Final remarks 

The construction of a microscopic field-theoretic framework to study the role of 
noise and dissipation in a realistic, yet efficient, way is still not accomplished. We 
believe this task will require the development of systematic (controlled) analytic 
approximations to simplify memory kernels, as well as effective numerical methods 
to deal with arbitrary colored noise. In the case of the chiral transition in heavy 
ion collisions, effects brought about by the expansion of the plasma and by its 
finite size will also bring corrections to this picture. Some of these issues will 
be addressed in future publications 122]. 
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